Solitons in the Salerno model with competing nonlinearities 
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We consider a lattice equation (Salerno model) combining onsite self-focusing and intersite self- 
defocusing cubic terms, which may describe a Bose-Einstein condensate of dipolar atoms trapped 
in a strong periodic potential. In the continuum approximation, the model gives rise to solitons in 
a finite band of frequencies, with sech-like solitons near one edge, and an exact peakon solution at 
the other. A similar family of solitons is found in the discrete system, including a peakon; beyond 
the peakon, the family continues in the form of cuspons. Stability of the lattice solitons is explored 
through computation of eigenvalues for small perturbations, and by direct simulations. A small 
part of the family is unstable (in that case, the discrete solitons transform into robust pulsonic 
excitations); both peakons and cuspons are stable. The Vakhitov-Kolokolov criterion precisely 
explains the stability of regular solitons and peakons, but does not apply to cuspons. In-phase and 
out-of-phase bound states of solitons are also constructed. They exchange their stability at a point 
where the bound solitons are peakons. Mobile solitons, composed of a moving core and background, 
exist up to a critical value of the strength of the self-defocusing intersite nonlinearity. Colliding 
solitons always merge into a single pulse. 

PACS numbers: 05.45.Yv, 63.20.Ry, 63.20.Pw 
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I. INTRODUCTION 

It is commonly known that the dynamics of nonlinear 
lattices are drastically different in generic nonintegrable 
systems, a paradigmatic example being the discrete non- 
linear Schrodinger equation (DNLS; see reviews 1]) and 
in exceptional integrable models, a famous example of the 
latter being the Ablowitz-Ladik (AL) equation [3. Only 
in the latter case are exact solutions for soliton families 
available. In nonintegrable systems solitons are sought in 
a numerical form 1] or by means of a variational approx- 
imation 3]. As the difference between the DNLS and AL 
equations is in the type of the nonlinear terms - onsite 
or intersite - it was quite natural to introduce a com- 
bined model that includes the cubic terms of both types, 
and thus allows one to consider a continuous transition 
between the AL and DNLS equations. The combined 
equation, known as the Salerno model (SM) 4], is 
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where is the complex field amplitude at the n-th site 
of the lattice, the over dot stands for the time derivative, 
and real coefficients ja and u account for the nonlinearities 
of the AL and DNLS types, respectively. 

The SM was studied in a number of other works, see 
Refs. and references therein. It is known that it 

conserves the norm and Hamiltonian, 
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In the above-mentioned works, it has been demon- 
strated that Eq. gives rise to static (and, sometimes, 
moving 0-0) solitons at all values of the DNLS param- 
eter u , and all positive values of the AL coefficient, ja. 
If u is negative, one can make it positive by means of 
the staggering transformation^ (— l)^^n, and then 

setting u = +1, by the rescaling ^n/V^ (unless 

u = 0). However, the sign of ja cannot be altered. In 
particular, the AL model proper {v = 0) with /i < does 
not give rise to solitons. The latter circumstance sug- 
gests considering soliton dynamics in the SM with /i < 0, 
i.e., with competing nonlinearities^ which is the subject 
of the present work. In this connection, it is necessary to 
stress that expressions (O and ^ for the SM's dynam- 
ical invariants remain valid if 1 + /i|^nP takes negative 
values at some sites, due to /i < 0. 

While the SM was originally introduced in a rather 
abstract context, it has recently found direct physi- 
cal realization, as an asymptotic form of the Gross- 
Pitaevskii equation describing a Bose-Einstein conden- 
sate of bosonic atoms with magnetic momentum trapped 
in a deep optical lattice |^]. In that case, the onsite 
nonlinearity is generated, as usual, by collisions between 
atoms, while the intersite nonlinear terms account for 
the long-range dipole-dipole interactions. Note that the 
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latter interaction may be attractive {fi > 0) or repul- 
sive (/i < 0), if the external magnetic field polarizes the 
atomic momentum along the lattice or perpendicular to 
it, respectively. 

This report is organized as follows. In Section ^ we 
develop a continuum approximation (CA), and investi- 
gate the corresponding solitons in an analytical form. It 
is found that, although they might exist in a semi-infinite 
band of frequencies, they actually occupy a finite band, 
with an exact peakon solution at the edge of the band. 
A family of discrete solitons is constructed by means 
of numerical-continuation methods in Section Hill They 
form a family of regular solitons, including a peakon, sim- 
ilar to what was found in the CA, but the lattice solitons 
extend beyond the peakon in the form of cuspons. In Sec- 
tion nil Al the soliton stability is explored by means of 
standard Floquet analysis (computation of perturbation 
eigenvalues or Floquet multipliers) and in direct simula- 
tions, with the conclusion that only a small part of the 
family is unstable. Two-soliton bound states are reported 
in Section HII Bl where it is demonstrated that stability 
exchange between in-phase and out-of-phase states oc- 
curs at a point where the bound solitons are peakons. 
Moving solitons are considered in Section IIVI where it 
is found that they exist up to a critical strength of the 
intersite self-defocusing nonlinearity, and collisions be- 
tween them always lead to fusion into a single soliton. 
The paper is concluded in Section IVl 



II. CONTINUUM LIMIT 



which may give rise to solitons, provided that a; < and 
< 1. (The absence of solitons for > 1 implies that 
if the intersite self-defocusing, accounted for by /i < 0, 
is stronger than the onsite self-focusing, the self-trapping 
of solitons is impossible in the CA). Equation Q can be 
cast in the form ((l/|/i|) - 1)"^ U^^ = -W'{U), where 
the effective potential is 
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the expansion of the potential ((Hj) for [/^ ^ is 1^ 
[-W^ + (1 - U^] /2 . This form of the equation 
shows that solitons exist in a finite band of frequen- 
cies, < < 1, rather than in the entire semi-infinite 
band, 1^ > 0, where the linearization of Eq. {TJ produces 
exponentially decaying solutions that could serve as the 
solitons' tails. The reduction of the semi-infinite band to 
a finite one is typical for soliton families in models with 
competing nonlinearities, such as the cubic-quintic NLS 
equation 10]. Further, it follows from the divergence of 
potential (|5|) at U'^ = l/|/i| that the solitons's amplitude 
A, which is a monotonously increasing function of 1^, is 
smaller than 1/ v4/? for < 1^ < 1, and A = 1/y4/? 

n = i. 

Solitons can be found in an explicit form near the 
edges of the existence band: at small (i.e., small Q), 

U{x) ~ Y^^^^4/T^^^^I/^^^c^ ^Y^2|a;|x^ , while precisely at 
the opposite edge of the band, 00 = 1 — l/|/i| (i.e., Q = 1), 
the exact solution is a peakon^ 



To introduce the CA in Eq. (pP), we define ^{x^t) = 
e^'^^(x,t), and expand ^n±i ~ ^±^a: + (l/2)^a:a:, where 
^ is now treated as a function of the continuous coor- 
dinate X, which coincides with n when it takes integer 
values. After that, the continuum counterpart of Eq. ^ 
is derived. 



i^t = -2{l-\^l\)\^f^-(l-\f^\\^f)^a 



(4) 



where we have set = +1 and /i < 0, as stated 
above. Equation ) conserves the norm and Hamilto- 
nian, which are straightforward counterparts of expres- 
sions and 
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Soliton solutions to Eq. Q are sought as \1/ = e {x), 
with a real function U obeying the equation 



(fU_ _ L0 + 2{l-\ii\)U^ 

1^ ~ 1 - 



(7) 



^peakon 



= (l/v4^) exp (-V(1/ImI)-1N) • (9) 



In other words, at a given frequency uo^ the peakon solu- 
tion is found at 



\^i\ = \^JiJ>\ = l/{l-Lo). 



(10) 



Note that norm ((SJ of the peakon is 7r^/[6Y^|/i|(l — |/i|)], 
and its energy is also finite. Close to this point, i.e., for 
< 1 — 1^ <C 1, the solution is different from the limiting 
form in a narrow interval |x| < (1 — 

where the peak is smoothed. 

Finally, the CA based on Eq. is valid if the intrinsic 
scale of all continuum solutions, that may be estimated 
through the curvatu re of the soliton's profile at x = 
as / ~ ^1 \/WxxW\^ is large, / ^ 1 (recall the lattice 
spacing is 1 in the present notation). According to Eq. 
©, the latter condition implies (l/|/i|) — 1 <C 1 (i.e., 
strictly speaking, the CA applies in the case when the 
competing nonlinearities in the SM nearly cancel each 
other). 

It is relevant to note that, in the usually considered 
"non-competitive" version of the SM, with /i > (i.e., 
self- focusing intersite nonlinearity), the CA yields soli- 
tons in the entire semi-infinite band, < 0. 
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III. STATIONARY DISCRETE SOLITONS 

In order to find discrete solitons in a numerical form, 
we looked for solutions to Eq. (^) which are localized 
and time periodic with frequency cj^ = 27r/T5 11] (that 
is related to uo in the continuum equation by a;^ = uo — 2). 
Soliton solutions of this form are widely known for the 
DNLS limit (/i = 0) and hence it is posible to make a 
numerical continuation of such solutions for /i < by 
adiabatic changes of the model parameter ji and succes- 
sive applications of the shooting method 0, For 
this purpose we consider soliton solutions of the above 
form as fixed points of the map Tt^^^ defined as 



(11) 



i.e Tt^^^ is the time-evolution operator over T5 of dynam- 
ics dictated by Eq. (^). Then, using usual techniques for 
finding fixed points of maps (such as e.g. the Newton- 
Raphson algorithm) we can find the numerically exact 
soliton for a given frequency 00^ and /i. In general all the 
soliton solutions were computed starting from the DNLS 
limit, /i = 0, and increasing at a fixed value of cob' The 
continuations were performed using the shooting method, 
with an increment S{\/i\) = 10~^ at each step, or smaller 
if higher accuracy was needed. 

As shown in the previous section, the soliton family 
in the continuum equation ends with the peakon 
solution 0. To compare the numerically determined 
shape of the discrete solitons with the feasible peakon 
limit, we fitted the solitons' tails to the asymptotic form, 
|^n| = ^exp(— r (|n — no|)), with constant A, F, and no, 
which follows from the linearized equation (^) for large 
|n|. This procedure yielded the decay rate, F = F(/i,cc;5), 
amplitude, A = A{/j.,ujij) (and the center's position no), 
as functions of parameters /i and 00^ of the soliton family. 
Once A{jii^u;i}) and no were found, we defined 7(/i,a;5) 
= A — \^no \ to measure a deviation of the true discrete 
soliton from a conjectured peakon shape obtained by for- 
mal extension of the tail inward. 

In Fig. ^a) we show the evolution of 7 produced by 
several continuations of the lattice-soliton solutions (at 
different frequencies coi)). We define jiipicob) as a value of 
/i at which an exact discrete peakon of internal frequency 
uji) is found, that we realize as vanishing of ^{fi^uoi)) at 
ji = jip. In Fig. n^b) we plot the evolution of the soli- 
ton's amplitude as the continuation is performed. It is 
observed that the amplitude increases with reaching 
the predicted value, l/y^^T^, at the exact peakon solu- 
tion. 

A noteworthy result, evident from Fig. ^ is the per- 
sistence of the discrete solitons beyond the peakon limit 
(which means continuability of the solutions to 7 < 0). 
The apparent intersection of different curves at one point 
in Fi^Ja) is a spurious feature (see the inset in the fig- 
ure): an accurate consideration shows that the curves 
actually intersect at close but different points. In con- 
trast, the intersection of the curves in Fig. ^b) indeed 




FIG. 1: (a) The mismatch with the peakon shape, 7 , as 
a function of for discrete solitons found at different fre- 
quencies ujb. (Note in the inset that there is no common in- 
tersection of all the curves), (b) The soliton's amplitude vs. 

The axes are rescaled to show that the amplitude of the 
peakon solutions (attained at \fi\ = |Mp|) are equal to l/^|/i|, 
as predicted by the continuum approximation. 



happens at a single point, which corresponds to the dis- 
crete solitons taking the peakon shape. 

Figure El displays typical examples of the numerically 
found discrete solitons. It demonstrates that the solu- 
tions corresponding to 7 < are cuspons^ with a super- 
exponential shape, that do not exist in the continuum 
equation The discrete character of the SM with the 
competing nonlinearities allows this new type of solution 
(as happens with the quasi-collapsing states in the stan- 
dard DNLS equation in two dimensions 0). Cuspon 
solutions continue into the region of > 1, where the 
CA yields no solitons, but, due to the sharp change of the 
solution with the increase of |//|, finding numerical solu- 
tions at larger values of becomes increasingly more 
difficult. 

In Fig. Ol^a) we compare the line of the existence of the 
peakons in the continuum limit, and the actual location 
of discrete peakons. It is seen that the agreement between 
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FIG. 2: Generic examples of three different types of discrete 
solitons, for Ub = —2.091: a quasi- continuous sech-like solu- 
tion at = 0.3, a peakon at = 0.956, and a cuspon at 
\fi\ = 2.64. 



the CA and numerical findings is good for smaller (in 
this case, the discrete solitons are broad), while at larger 
Icj^l the discrete solitons are narrow, hence the agreement 
with the CA deteriorates. 



A. Stability of the solitons 

We have performed stability analysis of the discrete 
solitons by computing eigenvalues (Floquet multipliers, 
Xp) for modes of small perturbations [13, within the 
framework of the linearized equation of the SM The 
integration of a basis of initial perturbations over a pe- 
riod T5 of the discrete soliton gives the Floquet matrix 
that in our case is real and symplectic so that its eigenval- 
ues come in quadruplets {A^, 1/Af, Af, 1/Af}. Then, 
the discrete soliton solution to Eq.Q is linearly stable 
if all the Floquet multipliers lie on the unit circle of the 
complex plane, |Ai?| = 1. It is found that the solitons are 
linearly stable along the whole continuation, except for 
a relatively small region, as shown in Fig. ^a). The en- 
tire instability island in the (|cc;5|, |/i|) plane is displayed 
in Fig. ini^b). Note, in particular, that the peakon and 
cuspon solutions are stable. The stability of the discrete 
solitons was also checked by direct simulations of per- 
turbed solitons, using the full equation which cor- 
roborated the predictions of the linear analysis. 

Direct simulations of the evolution of perturbed un- 
stable solitons, a typical example of which is displayed 
in Fig. ^b), show that, after a transient stage, a local- 
ized pulson (showing simultaneous width and amplitude 
oscillations) is formed. The pulsons are (quasi-)periodic 
in time, and persist indefinitely. This behavior resem- 
bles that found in the ordinary two-dimensional DNLS 
equation in quasi-collapsing states P^] . 



A necessary stability condition for soliton families 
in models of the NLS type may be provided by the 
Vakhitov-Kolokolov (VK) criterion if the norm Af of 
the soliton is known as a function of its frequency ujij, the 
solitons can be stable against small perturbations with 
real eigenvalues, provided that dJV/dcoi) < 0. Although 
the applicability of the VK criterion to the present model 
has not been proven (and counter-examples are known, 
when solitons predicted by the criterion to be unstable 
are actually stable ili]), it is relevant to test the criterion 
here, numerically computing Af {cOb) according to Eq. 
The result is that the VK criterion precisely explains the 
stability and instability of the discrete solitons, except 
for the cuspons (see below), as shown in Fig. ^c). 

A noteworthy feature of the M {uoh) dependence is a 
divergence of the total norm due to the infinite contri- 
bution of the central site to expression ^ in the case of 
the exact peakon solution, with |*^>nol^ = 
ample of the M {uo})) dependence showing the divergence 
is plotted in Fig. O As concerns the cuspons, whose 
amplitude exceeds the critical value, l/y^^, the norm 
(0) converges for them, and features a positive slope (see 
Fig. [51), dM I d(jj\) > 0. The VK criterion predicts in- 
stability in this case; however, the direct computation 
of the Floquet multipliers (see above), as well as direct 
simulations, reveal no instability of the cuspons. Thus, 
while the VK criterion is perfectly correct for regular soli- 
tons and peakons in the present model, it is irrelevant for 
cuspons, cf. the situation in Ref. [l^ . 



B. Bound states of solitons and their stability 

We have explored bound states of discrete solitons in 
Eq. (^). For this purpose, we performed numerical con- 
tinuation in /i, starting with the well known bound states 
of the standard DNLS equation, at /i = 0. In that limit, 
two different types of bound states are known, in-phase 
and TT-out-of-phase ones, which are represented, respec- 
tively, by even and odd solutions. It is well known that 
only the states of the latter type are stable (isj . 

The numerical continuation of soliton bound states was 
performed for pairs of identical discrete solitons of a given 
frequency cOb and different distances between them. The 
continuation of in- and out-of-phase bound states gives 
bound states of peakons, see Fig. IHlJa). The latter so- 
lution is found at exactly the same value, /i = /ip(cc;5), 
which gives rise to the single peakon. We have also ex- 
amined the stability eigenvalues for the computed solu- 
tions. A remarkable feature of the bound states observed 
with increase of is the stability interchange between 
the in-phase and out-of-phase states, as shown in Fig. 
IHI^b), which occurs precisely at /i = jipiuoi))^ regardless of 
the separation between the bound solitons. 
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FIG. 3: (Color online), (a) The value of at which the soliton assumes the peakon shape: the prediction of the continuum 
approximation, Eq. (fTn|) (solid curve), and numerical results for discrete solitons (dots). Also shown is a small region where the 
discrete solitons are found to be unstable (for that purpose, the vertical axis shows rather than |/ip|). The inset displays 
the relative difference between the numerically found values of and the prediction, provided by the continuum 

approximation, (b) Zoom of the area in the (Icj^I, |/ip|) plane where the instability island is located, (c) Norm of the discrete 
solitons vs. the frequency, for = 0.884. 



IV. MOVING SOLITONS 

The integrable AL model gives rise to both static and 
moving solitons. In the DNLS equati on, moving solitons 
l2nj and collisions between them ^ were also stud- 
ied, using various methods. Soliton motion was also a 
subject of analysis in the SM, which was usually treated, 
in this context, as a perturbed version of the AL equation 

As a part of the present work, we have performed nu- 
merical continuation of mobile discrete solitons from the 
DNLS limit, where, in turn, they were earlier obtained by 
means of a continuation procedure initiated with the AL 
solitons . The continuation procedure is the natural ex- 
tension of the one applied to static discrete solitons {i.e. 
based on computing fixed points of maps by means of the 
shooting method) and allows obtaining mobile solutions 
whose profiles are repeated (translated a number of sites) 
after an integer number of periods of the internal oscilla- 
tion. These states are composed of a traveling localized 
core and an extended background, = ^'^^^ + ^n^^^, 
see Fig. {7\ The background is a superposition of nonlin- 
ear plane waves, and its amplitude is related to the height 
of the corresponding Peierls-Nabarro (PN) barrier, which 
is defined as the energy difference between two static so- 
lutions for the soliton, with a fixed frequency cj^, one 
centered at a lattice site, with no = n, and the other at 
an intersite position, with no = n ± 1/2. 



The result obtained along these lines in the present 
model is that the mobile solitons can only be contin- 
ued up to a certain critical value, fi = /jjc{uJb)^ close to 
but smaller in absolute value than fip{uji})^ at which the 
static discrete soliton becomes a peakon. The Floquet 
stability analysis reveals that the extended background 
of the mobile solitons is subjected to modulational insta- 
bility. (However, this is too weak to manifest itself in the 
simulations and it is only noticeable by looking at the 
Floquet spectra when the amplitude of the background 
is very high). On the the other hand we do not observe 
any localized eigenvector with eigenvalue |Af| > 1 and 
thus the core is not affected by any unstable perturba- 
tion. The stability of mobile solutions is corroborated 
when simulations of the dynamics are performed allowing 
for interesting numerical experiments (see below). The 
background amplitude is a growing function of having 
a very sharp increase when approaches /i = jud^b)^ 
see FiglHJa). This behavior of the background amplitude 
suggests that the PN barrier also grows with and be- 
comes very high near the critical point. To check this ex- 
pectation, we have computed the height of the PN barrier 
for the same frequencies cOb for which the mobile solitons 
were numerically calculated, using the energy definition 
as in Eq. Figure [S^b) confirms that the PN barrier 
dramatically increases when the continuation approaches 
the critical point, /i = /ic(^6), although the PN barrier di- 
verges not exactly at this point, but rather at /i = /ip(a;5). 
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FIG. 4: (a) The absolute value of the Floquet multiplier, |Af|, 
for the linearization of perturbations around discrete solitons, 
is shown vs. , at several fixed values of the frequency (which 
are chosen so as to make the instability intervals well sepa- 
rated). The soliton is unstable if |Af| 7^ 1- (b) A robust 
pulson generated from an unstable soliton at — 0.922. 



where the soliton assumes the peakon shape. 

The strong dependence of the PN barrier on /i sug- 
gests a numerical experiment to test the behavior of mo- 
bile solitons when the lattice's pinning force suddenly 
changes. To this end, we took an initial mobile soli- 
ton at values of /i and a;^ for which the PN barrier is 
low. Then we monitored the evolution of the moving 
soliton following an instantaneous change in the nonlin- 
earity, /i ji ^ dji = ji' ^ which makes the PN barrier 
essentially higher than experienced by the original soli- 
ton. The numerical experiments are illustrated by Fig. 
ini We observe that the core of the mobile lattice soliton 
does not become pinned due to the increase of the PN 
barrier, but rather accommodates itself, with some radi- 
ation loss, into a broader state with a smaller amplitude, 
so that the PN barrier, as experienced by the new state 
for /i = /i , is low enough to allow the soliton to remain 
mobile. Besides that, we observe an increment in the 
core's velocity, so that the larger the jump of the PN- 
barrier's height the faster is the new moving state. The 
fact that the sudden increase of the PN barrier does not 
prevent the motion of the soliton reveals, on one hand, 
that the relation between PN barrier and mobility is far 
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FIG. 5: The norm of the discrete solitons vs. the frequency, 
for = 0.884. 



from trivial, and on the other hand, that mobility is quite 
a robust feature. 

Finally, we simulated collisions between identical lat- 
tice solitons moving in opposite directions. The results 
show that the colliding solitons always merge into a sin- 
gle localized state, which subsequently features intrin- 
sic pulsations. If the PN barrier is low, the emerging 
pulse can itself move in a chaotic way, due to interaction 
with the lattice phonon field (radiation) generated in the 
course of the collision. On the contrary, for values of ji 
and uoh at which the original solitons experience a high 
PN barrier, the finally generated single soliton is always 
strongly pinned. 

The most notable and generic feature of the collision 
manifests itself in the merger scenario. When the cores of 
the mobile solitons collide, sudden delocalization is first 
observed, with transfer of energy from the collision point 
to adjacent lattice sites. Then, almost all the energy 
is collected back at the collision spot, and thus a single 
localized state emerges. An example of the collision is 
shown in Fig. ^1 This scenario was observed in all 
simulations of the collisions. The appearance of pulsons 
as the product of soliton collisions, as well as the fact that 
they also appear as asymptotic states of the evolution 
of perturbed unstable solitons (see Section ITlI A|) . shows 
the ubiquity of this type of localized excitations in the 
present model. 



V. CONCLUSION 

We have introduced a combination of the DNLS and 
Ablowitz-Ladik models with competing nonlinearities 
(self-focusing onsite and self-defocusing intersite cubic 
terms). The model may describe a condensate of bosonic 
dipoles trapped in a strong lattice potential. First, it 
was shown that the continuum counterpart of the model 
gives rise to solitons in a finite band of frequencies, with 
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FIG. 6: (a) Profiles of typical in-phase (top) and out-of-phase (bottom) bound states of two peakons, with different distances 
between their centers, at Uh = —3.086 and |/i| = = 0.645 . (b) Absolute values of the Floquet multiplies that determine 
the stability of three bound states, with the same fixed frequency, Ub = —3.086, and different separations between the solitons. 
The in-phase (top) and out-of-phase (bottom) bound states are stabilized and destabilized, respectively, at the point where 
the bound solitons are peakons, see panel (a). Unstable states are less unstable (with smaller absolute values of the Floquet 
multipliers accounting for the instability) if the distance between the solitons is larger. 




FIG. 7: The real and imaginary parts of the lattice wave field 
in a moving discrete soliton, for cub = —2.24 and fi = —0.7. 



broad, small- amplitude NLS-like solitons at one edge of 
the band, and an exact solution in the form of a peakon 
at the other. Numerical analysis of the discrete sys- 
tem yields a similar family of solitons, which also in- 
cludes a peakon; however, the family continues beyond 
the peakon, in the form of discrete cuspons. Stability 
of the lattice solitons was investigated by means of the 
computation of their Floquet spectra, and in direct simu- 
lations. It was found that only a small part of the soliton 
family is unstable; the evolution of the unstable solitons 
leads asymptotically to pulsons, i.e. localized solutions 
where the width oscillates. In particular, peakons and 
cuspons are stable. Additionally, it was found that the 
Vakhitov-Kolokolov criterion precisely explains the sta- 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 



FIG. 8: The background amplitude (a) and the height of the 
Peierls-Nabarro barrier (b), as functions of |//|, for three mo- 
bile discrete solitons. 



bility and instability of regular solitons, but does not 
apply to cuspons, for which it erroneously predicts in- 
stability. Bound states of identical solitons were also in- 
vestigated, revealing a stability exchange: the in-phase 
and out-of-phase bound states, which are unstable and 
stable, respectively, in the DNLS limit, exchange their 
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FIG. 9: The time evolution of the real part of the lattice wave 
field in a mobile state. Initially, a moving soliton is created 
corresponding to = 0.8 and iUb = —2.24. Then, /i jumps 
instantaneously to = + Sfi. After the jump, followed by 
emission of some transient radiation, the moving core becomes 
broader and starts to move faster. In (a) = 0.84, and in 
(b) =0.89. 

stability character exactly at the point where the bound 
solitons are peakons. 

We have numerically computed mobile solitons with 
the same high precision as for the static ones. Their 
structure is that of an exponentially localized core on top 
of an extended background (which is a superposition of a 
finite number of extended plane waves). The amplitude 
of the background is correlated with the height of the 
computed Peierls-Nabarro barrier. The mobile solitons 
exist up to a critical value of the strength of the self- 
defocusing intersite nonlinearity, which is lower than (but 
close to) the corresponding value for which the peakon 



state exist. Collisions between solitons always lead to 

their merger into a pulson state. 
1 1 1 1 1 




time 

FIG. 10: Contour plots showing the evolution of the lattice 
field |$n| in three cases of collisions between identical solitons 
moving in opposite directions. The soliton's frequency is LUb — 
-2.11, and |//| = 0.6 (a), = 0.8 (b), and \fi\ = 0.9 (c). 



Acknowledgements 

The authors are pleased to acknowledge Franz G. 
Mertens for useful discussions. J.G.-G. appreciates hospi- 
tality of the Center for Nonlinear Studies at Los Alamos 
National Laboratory and acknowledges financial support 
from the MECyD through a FPU grant. B.A.M. ap- 
preciates hospitality of the Department of Condensed 
Matter Physics at the University of Zaragoza. Work at 
Los Alamos is supported by the USDoE. This work was 
partially supported by MCyT (Project No. BFM2002 
00113), DGA and BIFL 



9 



[1] M. Johansson, S. Aubry, Y. B. Gaididei, P. L. Chris- 
tiansen, and K. 0. Rasmussen, Physica D 119, 115 
(1998); P. G. Kevrekidis, K. 0. Rasmussen, and A. R. 
Bishop, Int. J. Mod. Phys. B 15, 2833 (2001). 

[2] M. J. Ablowitz and J. F. Ladik, J. Math. Phys. 17, 1011 
(1976). 

[3] B. A. Malomed and M. I. Weinstein, Phys. Lett. 220, 91 
(1996). 

[4] M. Salerno, Phys. Rev. A 46, 6856 (1992); see also En- 
cyclopedia of Nonlinear Science, p. 819 (ed. by A. Scott, 

New York, Routledge, 2005). 
[5] D. Cai, A. R. Bishop, and N. Gr0nbech- Jensen, Phys. 

Rev. E 53, 4131 (1996); K. 0. Rasmussen, D. Cai, A. R. 

Bishop, and N. Gr0nbech- Jensen, ibid. 55, 6151 (1997). 
[6] D. Cai, A. R. Bishop, and N. Gr0nbech-Jensen, Phys. 

Rev. E 56, 7246 (1997). 
[7] S. V. Dmitriev, P. G. Kevrekidis, B. A. Malomed, and 

D. J. Frantzeskakis, Phys. Rev. E 68, 056603 (2003). 
[8] J. Gomez- Gardehes, F. Falo, and L.M. Floria, Phys. Lett. 

A 332, 213 (2004); J. Gomez- Gardenes, L.M. Fiona, M. 

Peyrard, and A. R. Bishop, Chaos 14, 1130 (2004). 
[9] Z. D. Li, P. B. He, L. Li, J. Q. Liang, and W. M. Liu, 

Phys. Rev. A 71, 053611 (2005). 
[10] Kh. L Pushkarov, D. L Pushkarov, and L V. Tomov, Opt. 

Quant. Electr. 11, 471 (1979); S. Cowan, R. H. Enns, S. 

S. Rangnekar, and S. S. Sanghera, Canad. J. Phys. 64, 



311 (1986). 

[11] Discrete solitons are also widely referred to as discrete 
breathers because of the ocillatory dynamics of the real 
and imaginary part of the complex variables {^n}- 

[12] J.L. Marm, and S. Aubry, Nonlinearity 9, 1501 (1996). 

[13] T. Cretegny, and S. Aubry, Phys. Rev. B 55, R11929 
(1997). 

[14] E. W. Laedke , K. H. Spatschek, and S. K. Turitsyn, 

Phys. Rev. Lett. 73, 1055 (1994). 
[15] J.L. Marm, S. Aubry, and L.M. Fiona, Physica D 113, 

283 (1998). 

[16] M. G. Vakhitov and A. A. Kolokolov, Sov. J. Radiophys. 
Quant. Electr. 16, 783 (1973); see also a review: L. Berge, 
Phys. Rep. 303, 260 (1998). 

[17] L M. Merhasin, B.V. Gisin, R. Driben, and B.A. Mal- 
omed, Phys. Rev. E 71, 016613 (2005). 

[18] T. Kapitula, P. G. Kevrekidis, and B. A. Malomed, 
Phys. Rev. E 63, 036604 (2001); D. E. Pelinovsky, P. 
G. Kevrekidis and D. J. Frantzeskakis, Physica D 212, 1 
(2005). 

[19] M. J. Ablowitz, Z. H. Musslimani, and G. Biondini, Phys. 
Rev. E 65, 026602 (2002). 

[20] L E. Papacharalampous, P. G. Kevrekidis, B. A. Mal- 
omed, and D. J. Frantzeskakis, Phys. Rev. E 68, 046604 
(2003). 



